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Abstract 

For comparison of inclusive jet cross sections measured at hadron-hadron colliders to next-to- 
leading order (NLO) parton-level calculations, the energy deposited in the jet cone by spectator 
parton interactions must first be subtracted. The assumption made at the Tevatron is that the 
spectator parton interaction energy is similar to the ambient level measured in minimum bias 
events. In this paper, we test this assumption by measuring the ambient charged track momentum 
in events containing large transverse energy jets at y/s = 1800 GeV and y/s = 630 GeV and 
comparing this ambient momentum with that observed both in minimum bias events and with that 
predicted by two Monte Carlo models. Two cones in rj-<f) space are defined, at the same pseudo- 
rapidity, 7], as the jet with the highest transverse energy (Erp), and at ±90° in the azimuthal 
direction, (f>. The total charged track momentum inside each of the two cones is measured. The 
minimum momentum in the two cones is almost independent of E^p and is similar to the momentum 
observed in minimum bias events, whereas the maximum momentum increases roughly linearly 
with the jet E^p over most of the measured range. This study will help improve the precision of 
comparisons of jet cross section data and NLO perturbative QCD predictions. The distribution of 
the sum of the track momenta in the two cones is also examined for five different E^ bins. The 
HERWIG and PYTHIA Monte Carlos are reasonably successful in describing the data, but neither 
can describe completely all of the event properties. 

PACS numbers: 
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Jet production at hadron colliders, the highest energy probe in particle physics, has been 
used to measure parton distribution functions, the running of the strong coupling constant, 
a s , and to search for new physics. At the Fermilab pp collider, the jet production rate 
has been measured for jets of 15-450 GeV at ^ = 1800 GeV Q, Q, Q, Q and jets of 
15-150 GeV at y^i = 630 GeV The production of jets involves the interaction of 

an individual parton (quark or gluon) from one beam hadron with a parton from the other 
beam hadron. Each of the interacting partons carries only a fraction of the parent hadron's 
momentum with the residual momentum remaining with the other (spectator) constituents 
of the hadron. In addition, there are interactions between the spectator constituents of the 
two hadrons which normally occur at low momentum transfers. Measurements involving the 
observed jets are compared to perturbative QCD predictions. For NLO perturbative QCD 
predictions, only the parton level cross section, i.e. the cross section of two partons producing 
either two or three partons in the final state, is calculated. After convolution with the parton 
distribution functions, this cross section is directly compared with experimental data. For 
these comparisons to be valid, the energy from spectator interactions, which may fall in the 
jet cone, must be subtracted from the experimentally observed jets. In hard interaction jet 
events, the energy outside the two primary jets consists of energy from spectator interactions 
(soft and semi-hard), initial and final state radiation and any hadronization leakage from the 
jet cones. Initial and final state radiation effects are part of higher order perturbative QCD 
calculations and at least a portion of these effects are already included in NLO calculations. 

Because the CDF detector measures the momenta of low tracks more accurately than 
the calorimeter measures their energies, we choose to work with the track momenta in our 
analysis. We will call the charged track momenta associated with spectator interactions the 
underlying event momentum. It is the momentum in an event which is not directly related 
to the hard interaction. Clearly, this is a working definition as a coupling exists between 
all aspects of a pp interaction. For example, the hadronization of the partons from the 
hard interaction and from the spectator interactions are ultimately linked as the final state 
hadrons must be colorless. In current QCD studies at hadron colliders, the underlying event 
energy in the jet events is assumed to be well approximated by the ambient energy in the 
events collected with minimal trigger requirements. Normally, these minimum bias events 
are triggered by presence of particles away from the beam in the forward and backward 
direction. The subtraction of the underlying event energy leads to the largest uncertainty 
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in jet cross section measurement for < 50 GeV |l[. A precise measurement of the 
spectator interaction energy is essential for the modeling/ understanding of non-perturbative 
QCD effects and for any quantitative improvement of the jet studies. Another important 
question is whether the presence of a hard interaction in the event influences the spectator 
interactions. 

The measurement of the momentum in minimum bias events is important in its own right 
as it is used to estimate the effect of pile-up events on any signal at hadron colliders, where, 
due to high instantaneous luminosity, several interactions may occur in the same bunch 
crossing. In this paper, we present a measurement of the momentum deposited far from the 
jets in pp interactions at ^/s = 1800 and y/s = 630 GeV and compare our measurement 
with the momentum observed in minimum bias events and with the predictions from two 
Monte Carlo models. The jet samples used in this analysis are the same as for the inclusive 
jet cross section measurements at the two center-of-mass energies. The study of interjet soft 
gluon radiation is also of special interest in QCD as its emission originates from the flow of 
color between jets. The analysis of such observables may lead to a better understanding of 
color neutralization 7J, r|. 

Q 

The study reported in this paper is complementary to our previous analysis [9J], which 
examined the evolution of event structure in low to moderate Et events in pp interactions at 
\fs = 1800 GeV by studying charged particle jets from 0.5 GeV/c to 50 GeV/ c. The previous 
study found that the momentum transverse to the leading jet rises rapidly in the 0.5-5.0 
GeV/c range and is almost constant when the leading charged particle jet has transverse 
momentum greater than about 10 GeV/c. 

To study the underlying momentum in jet events, we define two cones with radius R = 
y/(Ar]) 2 + (A0) 2 = 0.7 centered at 

Y] = T j( 1 )^ a nd (p = ( p( 1 ) ±90° where (t/ 1 -', 4 1 ) * s the centroid 
of the highest energy jet in the event as shown in Figure ^ The sum of the transverse 

, r ii , i • ji , i l i j 7->90, rain j j-)90,max i r->90,max ■ 

momenta of all tracks m the two cones is labeled F T and f T , where f T is 
higher of the two values. By definition, p®°' max should contain a larger contribution from 
initial and final state radiation than p^°> 7mn , j n ^ e approximation of a negligible four parton 
final state component, p^°' mm [ s a measure of the underlying momentum in the jet event. 
In minimum bias data, we perform a similar analysis but with the cone centroid selected 
randomly in the central rapidity region, \r]\ < 0.5. We also use a second procedure, the Swiss 
cheese method, in which the transverse momenta of all the tracks except those in the two 
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90 cones 
used in this 
analysis 




Lead Jet 



Second Jet 



FIG. 1: An example of a two jet event in the detector region under study. The cones used for 
the determination of the underlying event contribution are at r/ = 77' ^ and 4> = (f>W ± 90° where 
(r/ 1 ), 4>^) is the centroid of the highest Et jet in the event. 



or three highest energy jets are summed and compared with Monte Carlo predictions and 
minimum bias data. A study of this type was first suggested in Ref. jlfll |. 

The data were collected using the CDF detector [ll| with the Fermilab Tevatron Collider 
at y/1 = 1800 GeV (1994-1995) and = 630 GeV (1995). The CDF detector is a multi- 
purpose detector consisting of a tracking system in a solenoidal magnetic field, calorimeters, 
muon chambers and two arrays of scintillator counters (BBC) located at ±5. 8m from the 
nominal interaction point along the beam direction, covering the 3.2 < \rj\ < 5.9 region. 
Minimum bias events were triggered by a coincidence of hits in these counters. The BBC 
cross section is 51.15±1.60 mb compared to a total inelastic cross section for pp interactions 
of 60.33 ± 1.40 mb at ^/s = 1800 GeV l^J. The jet data were collected using four trig- 
gers requiring a cluster of energy in the calorimeter with Et > 20, 50, 70 and 100 GeV at 
y/s = 1800 GeV. These data samples, and the jet clustering and energy corrections have 
been described in detail in |1] . The data at y/s = 630 GeV employed two triggers requiring 
a cluster with Et > 5 and 15 GeV respectively. The jet energies were corrected for any 
energy loss in the detector. At both energies we use only those events in which the centroid 
of the highest energy jet is within the central rapidity region, \r]\ < 0.5. 

The tracking system consists of a silicon vertex detector (SVX'), a vertex tracking cham- 
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ber (VTX), and a central tracking chamber (CTC) Q]. The vertex reconstruction is 
performed using information from the VTX and the CTC. In this analysis, the jet events 
were required to have one and only one primary vertex of high quality (corresponding to 
a high track multiplicity). For inclusive jet analyses in general, there is no restriction on 
the number of vertices as long as there is at least one high quality vertex. The one vertex 
requirement in this analysis is implemented in order to restrict the events to those in which 
only one interaction occurred during that beam crossing. For minimum bias data, the re- 
quirement is changed to one vertex of medium quality (corresponding to a lower minimum 
track multiplicity, but one resulting from beam-beam rather than beam-gas interactions). 

Track reconstruction takes place primarily using hit information from the CTC. In order 
to ensure a high quality for the reconstructed tracks, each track is required to have at least 
four hits in each of the five axial super-layers and hits in at least one stereo super-layer. 
The momentum resolution in the rapidity region \rj\ < 1 is better than SPt/Pt 2 < 0.002 
(GeV/c)" 1 . We require the tracks to have Pt > 0.4 GeV/c and to be within 5 cm in the 
longitudinal and 0.5 cm in the transverse direction of the pp vertex. The uncertainty on the 
quantities measured in this analysis is evaluated by loosening these cuts to 10 cm and 5 cm 
respectively. The charged track reconstruction efficiency is uniform in rapidity for \rj\ < 1 
and is on average ~ 92 ± 3% Q. Th e efficiency drops to 80% at low P T (0.4-0.5 GeV/c) 
and to 60% in the region 1.0 < \q\ < 1.2. We correct the data for the inefficiency in both 
regions. The main systematic uncertainty in our analysis arises from the track selection 
criteria and from the track reconstruction efficiency. The data were compared to Monte 
Carlo predictions from the programs HERWIG Q and PYTHIA Q- At yfs = 1800 GeV, 
four samples of jet events were generated with HERWIG and PYTHIA, with a minimum 
transverse momentum for the hard scattering of 20, 40, 60 and 80 GeV, for the four samples. 
The leading jet in the generated distributions was required to have a transverse energy of 40, 
75, 100 and 130 GeV, respectively. The output from both Monte Carlo programs consists 
of the 4-vectors of the final state hadrons. For comparison to the Monte Carlo predictions, 
the data were corrected for the track reconstruction efficiency. 

In Figure El a comparison of the p® ' max / mm distributions between data and Monte Carlo 
is shown. p®°' max increases as the leading jet Et increases, both in the data and in the Monte 
Carlo predictions. The p^°' mm distributions are almost independent of E^p indicating that 
any contribution from higher order radiation, at least in this r]-(f) region, is small. Both 
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FIG. 2: P; 
jet at yfi = 1800 GeV. 



and their difference APjP as a function of the Et of the highest energy 



the Monte Carlo predictions and the data show a similar behavior. The average values of 
p90,max p90,mm given in Table for different intervals of Elp , for the data and 
for the two Monte Carlos. Good agreement is observed with HERWIG, while PYTHIA 
lies above the data. The parameters of the underlying event model in PYTHIA can be 
adjusted more easily than those in HERWIG. We have attempted to reach a better agreement 
with the predictions from PYTHIA by using a more modern parton distribution function 
(CTEQ4L instead of MRSG Q), 

using the option of varying impact parameters with a 
matter distribution inside the hadron described by a simple Gaussian (MSTP(82)=3), and 
by decreasing the regularization scale of the transverse momentum spectrum for multiple 
interactions (Pto) to 2.0 GeV/c from the default value of 2.3 GeV/c. (Such a decrease causes 
the double parton scattering component of the underlying event to be less hard, leading to a 
better agreement with the data.) Table ITT1 summarizes the Monte Carlo's default and tuned 
parameters. The behavior of PYTHIA with the adjusted parameters can be observed in 
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TABLE I: Average Pt inside the max and rain cone at rj = rp V and (f) = 4> {l) ± 90° for yfl = 1800 
GeV data. In data, the first errors shown are statistical and the second are systematic. 



4 1 ' 


DATA 


HERWIG 


PYTHIA 


PYTHIA 


(GeV) 






(default) 


(tuned) 


r>90, max t t / \ 

P T (GeV/c) 


40-80 


2.04±0.09±0.21 


1.92±0.04 


2.43±0.04 


2.19±0.04 


80-120 


2.64±0.09±0.19 


2.49±0.05 


3.39±0.06 


2.96±0.06 


120-160 


2.89±0.09±0.22 


2.95±0.05 


3.69±0.06 


3.56±0.06 


160-200 


3.27±0.10±0.22 


3.21±0.07 


4.02±0.09 


3.93±0.09 


200-270 


3.64±0.21±0.24 


3.59±0.16 


4.35±0.17 


4.24±0.19 


p 90,mm (GeV/c) 


40-80 


0.37±0.03±0.06 


0.38±0.01 


0.54±0.01 


0.38±0.01 


80-120 


0.47±0.02±0.07 


0.43±0.01 


0.61±0.01 


0.44±0.01 


120-160 


0.42±0.02±0.06 


0.45±0.01 


0.64±0.01 


0.48±0.01 


160-200 


0.46±0.02±0.06 


0.48±0.01 


0.62±0.02 


0.53±0.02 


200-270 


0.53±0.05±0.07 


0.50±0.02 


0.63±0.03 


0.53±0.04 



Figure El and in the last column in Table |U The tuning leads to a better agreement with 
the data in the low E^p region but leaves PYTHIA still somewhat larger in the high E^ 
region. 

In Figure EJ the total charged track momentum in the two cones (min+max) is shown 
for five different bins of E^. The effects due to large angle (away from any jet) soft gluon 
emission are expected to be appreciable when the transverse momentum in the min+max 
cones (p"!p ax + p™ m ) is larger than a few GeV and when the ratio of the lead jet transverse 
momentum to the transverse momentum in the cones is large. Such emissions are included 
in an approximate way in existing parton shower Monte Carlos and a detailed comparison 
may lead to improvements in their treatment 

Q- Qualitatively, HERWIG and PYTHIA 
agree with the data, although the PYTHIA prediction tends to be slightly higher than the 
data for larger values of transverse momentum in the two cones |2l| . 

The average number of tracks found inside the two cones is shown in Figure plotted as 
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TABLE II: Monte Carlo's default and tuned parameters. MSTP(82) defines the structure of the 
multiple parton interactions; PARP(82) is the regularization scale of the transverse momentum 
spectrum for multiple interactions (with MSTP(82)>2); PARP(85) and PARP(86) are the prob- 
ability that the multiple interaction produces two gluons with color connections to the nearest 



neighbors (or as a closed gluon loop) 





PYTHIA 


HERWIG 




default 6.115 


tuned 6.115 


tuned 6.115 


default 5.6 




1800 and 630 GeV 


1800 GeV 


630 GeV 


1800 and 630 GeV 


Parton distribution function 


MRSG 


CTEQ4L 


CTEQ4L 


CTEQ3L 


MSTP(82) 


1 


3 


3 




PARP(82) 




2.0 GeV/c 


1.4 GeV/c 




PARP(85) 


0.33 


1 


1 




PARP(86) 


0.66 


1 


1 





a function of E T . A slightly higher track multiplicity is observed in both of the simulations 
compared to the data. 

In Table nm the mean values of the total track Pt and the mean number of tracks inside 
a cone randomly placed in the region \rj\ < 0.5 are shown for all minimum bias events and for 
those with a high quality vertex only. For the entire sample, the mean transverse momentum 
(P^ fB,cone ) in the cone is about 0.36 ± 0.04 GeV/c, while restricting the sample to events 
having a high quality vertex, the transverse momentum increases to 0.57 ±0.06 GeV/c. The 
average for the p^°' min cone over the measured Erp range is approximately 0.45 GeV/c, or 
midway between the above values. 

The track multiplicity and track momentum distributions in minimum bias data are 
shown in Figure El with the number of entries in the simulation normalized to the number 
in the data. The transverse momentum distribution at high Pt is not well-reproduced by 
HERWIG, which has virtually no tracks with Pt > 4 GeV/c. The absence of high Pt tracks 
indicates the lack of a semi-hard processes in the HERWIG model of minimum bias events. 
In contrast, PYTHIA reproduces the transverse momentum distribution considerably bet- 
ter. The model of multiple parton interactions incorporated in the PYTHIA description of 
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minimum bias events |19| allows for the possibility of high transverse momentum tracks. 
Neither HERWIG nor PYTHIA appears to correctly describe the high multiplicity end of 
the track multiplicity distribution. 

As previously described, the Swiss cheese distribution is formed by summing the trans- 
verse momentum of the tracks in the central region < 1), excluding the transverse 
momentum of the tracks in a radius 0.7 from the center of the two (or three) most energetic 
jets in the event (where an requirement of 5 GeV has been placed on each jet). 

The sum of the track transverse momentum in the central region for the 2-jet subtracted 
and 3-jet subtracted distributions is shown in Figure[7|for the data, HERWIG and the version 
of PYTHIA tuned for best agreement with the max/min cone data. In the simple picture 
presented earlier, the difference between the Swiss cheese level with the 2 highest Et jets 
subtracted and the corresponding minimum bias level should be proportional to the NLO 
(third parton) and higher order contributions. The Swiss cheese level with the three highest 



• MAX Data 




•> t !■ MAX-MIN Data 
X » MAX-MIN HERWIG 



■ MAX-MIN PYTHIA6.1 1 5 (tuned) 



5 
< 




1 25 50 75 100 125 150 175 200 225 250 

E T of leading jet (GeV) 



FIG. 3: P* ) ' maa! J p^ mm and their difference AP T ° function of the Et of the highest energy 
jet at yfs = 1800 GeV. PYTHIA has been tuned to reproduce the data. 
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FIG. 4: The distributions for the total pt in the sum of the max and min cones is plotted for five 
different bins of the Er of the highest energy jet. Data, HERWIG and PYTHIA distributions are 
shown at y/s = 1800 GeV. 

Et jets subtracted should have little or no NLO contribution. Both the 2-jet subtracted and 
3-jet subtracted distributions increase as the lead jet E T increases, with the slope being less 
for the 3-jet subtracted case. The 3-jet subtracted Swiss cheese average Py/(unit rj— (f>) is 
0.92±0.09 GeV/c compared to 0.37±0.04 GeV/c observed in minimum bias data with high 
quality vertices and 0.23 ± 0.04 GeV/c in all minimum bias events. The larger momentum 
observed in the 3-jet subtracted Swiss cheese distribution indicates additional contributions 
than just those from the soft underlying event. These contributions include hadronization 



221, as 



from the jets (splash-out), double-parton scattering, higher order radiation effects 
well as contributions from 3 rd jets that fail the Et threshold cut of 5 GeV. For comparison, 
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FIG. 5: Number of tracks in the max and min cone as a function of the Et of the leading jet. 
Data, HERWIG and PYTHIA distributions are plotted at yfs = 1800 GeV. 

the average momentum/ (unit r)-<j)) in the min(max) cone is 0.29±0.04 (1.91±0.14) GeV/c. 

To study the energy dependence, we have analyzed jet and minimum bias data at ^/s = 
630 GeV. In Figure El p®°' max and p^°> mm (and their difference) are plotted as a function 
of The y/s = 630 GeV data shows a similar behavior as was observed at = 1800 

GeV but the overall magnitudes are lower. The average p^°' mm (Py ' maa: ) ^ ^ _ q^q 
GeV is 0.25 ± 0.04 (1.43 ± 0.12) GeV/c, ~ 0.2 (~ 1.5)GeV/c lower than what is observed 
at v/i = 1800 GeV. Both HERWIG and PYTHIA reproduce the data at 630 GeV well. 
PYTHIA has been tuned as for the analysis at 1800 GeV, but with the regularization 
scale, Pro, set to 1.4 GeV/c. A dependence of Pro on the center of mass energy has been 
implemented in versions of PYTHIA after 6.12 according to the model described in j^. 
This model, however, predicts a value for P^o of 1.9 GeV/c for the 630 GeV data while 
in the 1800 GeV data the prediction is 2.3 GeV/c, showing a smaller dependence on the 
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FIG. 6: Distribution of track multiplicity (left) and transverse momentum (right) in the \fs = 1800 
GeV minimum bias sample. The inset in the right plot shows the low Pt portion of the spectrum 
on a linear scale. 

center-of-mass energy than that observed in our data. Using the default values the PYTHIA 
predictions underestimates the number of charged particles at y/s = 630 GeV. 

The Swiss cheese distributions at 630 GeV are shown in Figure El A very good agree- 
ment between data and both Monte Carlos is observed if the three most energetic jets are 



TABLE III: Mean p* 1D > cone an d the mean number of tracks in a random cone of radius 0.7 in 
y/s = 1800 GeV minimum bias data. Only systematic errors are shown. Statistical errors are less 
than 0.5%. 







pMB,cone 
T 


Track 






(GeV/c) 


Multiplicity 


DATA 


all vertices 


0.36 ± 0.04 


0.45 ± 0.06 




high quality vertex 


0.57 ± 0.06 


0.69 ± 0.09 


HERWIG 


0.31 


0.44 


PYTHIA (tuned) 


0.35 


0.44 
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FIG. 7: P T um (Swiss cheese). The two and three most energetic jets in each event are subtracted 
from the total transverse momentum in the central detector region. Data, HERWIG and PYTHIA 
results are shown at y/s = 1800 GeV. 



subtracted. In the case where the two most energetic jets are subtracted, both Monte Carlos 
lie above the data, as was also observed at 1800 GeV. Again, the momentum is larger when 
three jets are subtracted {Pt/ '(unit rj— (f>) is 0.52±0.05 GeV/c ) than in minimum bias events 
with a high quality vertex(P^/ (unit rj-<p) is 0.34±0.03 MeV/c). 

In Table IIVI are shown the average value of the total track Pt and the mean number 
of tracks inside a cone in the central rapidity region in minimum bias data at \fs = 630 
GeV. The sum of the track transverse momenta is 20% lower with respect to the 1800 GeV 
sample. Figure ITUl shows the track multiplicity and momentum distributions for minimum 
bias events. Again, the number of entries in the simulation is normalized to the number of 
entries in the data. The track multiplicity distribution and the mean P T B,cone y dominated 
by the low edge of the steeply falling spectrum, is well reproduced by both Monte Carlo 
generators. Unlike the situation at 1800 GeV, PYTHIA fails to produce enough high Pt 
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FIG. 8: p9°>™ n ) P ^,max &nd theiv difference Ap 90 at ^ = 630 GeV as a function of the E T of 
the leading jet. 

tracks, although it still produces considerably more than HERWIG. 

TABLE IV: Data and simulation comparisons for minimum bias events at ^fs = 630 GeV. Average 
Prp ' and the average number of tracks in a random cone of radius 0.7 are shown. Only 
systematic errors are shown. Statistical errors are less than 0.2%. 





pM B ,cone 
T 

(GeV/c) 


Track 
Multiplicity 


DATA 


all vertices 
high quality vertex 


0.29 ± 0.03 
0.52 ± 0.05 


0.37 ± 0.05 
0.65 ± 0.08 


HERWIG 


0.26 


0.36 


PYTHIA 


0.28 


0.38 



In summary we have studied the momentum deposited in two cones at ±90° to the 
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FIG. 9: p^ um (Swiss cheese). The two and three most energetic jets in the events are subtracted 
from the total transverse momentum in the central detector region. Data, HERWIG and PYTHIA 
results are shown at yfs = 630 GeV. 

highest E T jet in hard interaction events at y/s = 1800 and 630 GeV. The maximum of 
the two cone energies increases with highest Et jet in the event whereas the minimum is 
flat. Both HERWIG and PYTHIA exhibit the same behavior but HERWIG provides a 
better description of the CDF data. The momentum in the min cone is midway between 
the levels observed in generic minimum bias events and minimum bias events selected by 
high track multiplicity. In the HERWIG minimum bias model, the generated tracks are 
too soft, and semi-hard or hard interactions should be added to the minimum bias events 
in order to better reproduce the data. PYTHIA, however, with an adequate tuning of 
its parameters, reproduces the charged particle distribution better than HERWIG for the 
1800 GeV minimum bias data, but is less successful at 630 GeV. We have measured the 
y/s dependence of the underlying event momentum in jet events and ambient energy in 
minimum bias events. The underlying momentum at y/s = 630 GeV is 45% lower than what 
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FIG. 10: Distribution of track multiplicity (left) and transverse momentum (right) in the -^s = 630 
GeV minimum bias sample. The inset in the right plot shows the low Pt portion of the spectrum 
on a linear scale. 

is observed in 1800 GeV data. These measurements will allow for more precise tunings of 
both the underlying event in Monte Carlo programs and the mechanisms for gluon radiation, 
help to reduce the uncertainties in future jet studies at the Tevatron and will lead to a better 
prediction of physics signals and backgrounds at the Large Hadron Collider. 
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